Search for H dibaryon on the lattice 
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We investigate the H-dibaryon, an I{J^) = O(O^) with s — —2, in the chiral and continuum 
regimes on anisotropic lattices in quenched QCD. Simulations are performed on modest lattices 
with refined techniques to obtain results with high accuracy over a spatial lattice spacing in the 
range of Us ~ 0.19 — 0.40 fm. We present results for the energy difference between the ground state 
energy of the hexa-quark stranglet and the free two-baryon state from our ensembles. A negative 
energy shift observed in the chirally extrapolated results leads to the conclusion that the measured 
hexa-quark state is bound. This is further confirmed by the attractive interaction in the continuum 
limit with the observed H-dibaryon bound by 47 ± 37MeV. 

PACS numbers: ll.15.Ha, 14.40.Lb, 12.38.Gc, 12.38.Mh 



I. INTRODUCTION 



Search for dibaryons is one of the most challenging 
theoretical and experimental problems in the physics of 
strong interaction. In the non-strange sector, only one 
dibaryon, the deuteron, is known experimentally. In 
the strange sector, on the other hand, it is still un- 
clear whether there are bound dibaryons or dibaryon 
resonances. Among others, the flavor-singlet state (uud- 
dss), the H-dibaryon, has been suggested to be the most 
promising candidate [3. The H-dibaryon may also be a 
doorway to strange matter that could exist in the core of 
neutron stars and to exotic hyper- nuclei 0, Q • A deeply 
bound H with the binding energy Bh > 7 MeV from the 
AA threshold has been ruled out by the discovery of the 
double A nuclei [1] . The double- hypernucleus events have 
either more than one interpretation for the products or 
the possibility of production of excited states. The anal- 
ysis of the events in [3| ignored the possibility that the 
single or double hypernucleus was produced in an ex- 
cited state in which case the value of the binding energy 
would increase by the excitation energy. Thus there still 
remains a possibility of a shallow bound state or a reso- 
nance in this channel. Since the Bag Model prediction of 
a deeply bound H-dibaryon, with large binding energy of 
0(100) MeV 0, many experiments have been triggered 
to look for this possible particle, but few of them have 
confirmed the existence of the H-dibaryon [6l-[l0j . The ex- 
periments were inspired by the Skyrme model prediction, 
and the experimental discoveries have in turn spawned in- 
tense interest on the theoretical side, with studies ranging 
from chiral soliton and large Nc models, quark models. 



phase-shift analysis and QCD sum rules jllj . Summa- 
rizing the previous theoretical and experimental investi- 
gations a slightly bound or unbound H-dibaryon is pre- 
dicted. 

In the search for such exotic states lattice QCD 
plays an important role in which precise predictions for 
hadronic observables with quantifiable uncertainties are 



made. A considerable interest in H-dibaryon started with 
Jaffe's work [Bj demonstrating the role of the chromo- 
magnetic interaction in the stability of light multiquarks. 
Since then a number of quenched LQCD calculations 
have been performed for the search of the H-dibaryon 
but no definite conclusions have been reported. Earlier 
lattice investigations [l^, gave somewhat mixed and 
contradicting results on the spatially localized resonance 
status of H-dibaryon. These studies however, suffered 
from low statistics, relatively large quark masses and con- 
siderable finite-size effects could not be rules out for the 
smaller lattice size. More precise studies on large vol- 
umes concluded an unbound H-dibaryon in infinite vol- 
ume limit while others reported hints of a bound 
H-dibaryon for a range of light-quark masses [l3| • 

Very recently, NPLQCD and HALQCD Cohaborations 
reported results from their fully dynamical lattice calcu- 
lation that shed new light on the status of H-dibaryon 
[H, NPLQCD CoUaboration presented a strong 

evidence for a bound H-dibaryon from their calcula- 
tions performed in four lattice volumes. Using Lusher 
method [10] to extract two-particle scattering ampli- 
tude below threshold, NPLQCD collaboration found H- 
dibaryon bound by 16.6 MeV at a pion mass of m^r = 390 
in the infinite volume limit [l^. Using their recently 
proposed approach of baryon-baryon potential [2l| , the 
HALQCD Collaboration performed calculations in three 
lattice volumes at three different quark masses and re- 
ported a bound H-dibaryon with the binding energy of 
30 - 40 Mev for the pion mass of 673 - 1015 MeV |ll|. 
These calculations provide strong evidence of capability 
of Lattice QCD in calculating the energy of simple nuclei, 
with H-dibaryon being an example. Having said that, it 
still remains an open question whether the H-dibaryon is 
bound at the physical point and with the inclusion of the 
electroweak interactions. This provides a strong motiva- 
tion for pursuing numerical calculations at smaller lat- 
tice spacings, and over a range of quark masses including 
those of nature. 

The target of this work is to address the status of the 
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H-dibaryon by calculating the mass differences between 
the candidate H-dibaryon and the free two-baryon states 
in the continuum limit and at physical quark masses. 
Using some refined methods and techniques, we carry a 
multi-lattice spacing analysis at and near physical pion 
mass on improved anisotropic lattices with an attractive 
feature that with modest lattice sizes one can access large 
spatial volumes while having a fine temporal resolution. 
Rather than extracting the hadron mass from the ra- 
tio of two temporal nearest correlators, the Levenberg- 
Marquardt algorithm is adopted to solve the hyperbolic- 
cosine ansatz of hadron correlation functions. This is 
very useful in finding a larger temporal fit range, hence 
more clear signals for precise hadron masses. Continuum 
limit is also considered in this work which will provides 
a real physical picture of the H-dibaryon. 

The rest of the paper is organized as follows. The 
technical details of the lattice simulations are discussed 
in Sec. H, where we outline the construction of the cor- 
relation functions from interpolating operators and the 
actions used in this study. This section also discusses the 
procedures of chiral and continuum extrapolations. The 
results are presented and discussed in Sec. HI, where 
we attempt to take the chiral and continuum limits and 
address chiral, finite-spacing and quenching effects. This 
sets the stage for a discussion of lattice resonance signa- 
ture of the H-dibaryon lying lower than two-A channel 
masses in the physical regime. Finally, we present our 
conclusions in Sec. IV. 



II. SIMULATION DETAILS 

A. Choice of Interpolating fields 

The explicit construction of the operator for H- 
dibaryon requires the symmetrization of the colour and 
spinor indices of two triplets of quarks in order to obtain 
colour and spin singlet. Our choice for the appropriate 
operator is motivated by possible structure of H-dibaryon 
and based on the idea of diquark formulation and has the 
following form [H, [2^ : 



Oh{x) 
(udsuds) 



^{udsuds) — 3{ussudd) — 3{dssduu) 



(1) 



where the roman letters denote the colour indices, greek 
letters represent the spinor indices and e"*"^ the usual an- 
tisymmetric tensor defined over the range of their indices. 
Taking the symmetry properties of the e-tensor and the 
(C75)-matrix under the interchange of two indices into 
account, the H-dibaryon correlation function can be ob- 
tained from 



C„{x,t) =<OH{x,t)oU0)> 



(2) 



The hadron masses Mh and Ma, needed to obtain the 
energy shift AE = Eh — 2m\, are calculated by the fitting 



the correlation functions with a multi-hyperbolic-cosine 
ansatz 



C{t) ^Y,A,cosh[m,{T/2 - t)], 



(3) 



i=0 



where is the effective mass of the ith excited state 
and Ai the amplitude corresponding to this state. The 
method of calculation is straightforward in principle, 
not differing essentially from the calculation of hadron 
masses. We determine the mass of the ground state 
for each particle, and the mass difference AM{H — 2A) 
from the fit. In order to reduce the contaminations of 
the excited states, the maximum time separation is used 
to extract the results. This is achieved by adopting 
the Levenberg-Marquardt algorithm to solve this non- 
linear least-squares fit. The hyperbolic-cosine fits are 
performed over the time interval in which an acceptable 
value of the probability, used to estimate the goodness-of- 
fit of the data, is obtained. Considering the contribution 
of ground state only, the correlation function is fitted by 
the form 



C{t) = ^0 cosh[TOo (r/2 - i )] . 



(4) 



To account for the strong correlation of data in time, 
we use the full covariance matrix to construct the 
function 

= 5^[C(t„^o,mo)-Z?(tO]Mi7'[Cfe:^o,mo)-i?(tj)] 

(5) 

and obtain the the covariance matrix My as 



M,. 



1 ^= 



(6) 

where the D{k,ti) is the fcth correlator and D{ti) the 
mean value of the correlator at time ti, Nc denotes the 
total number of configurations. At minimum x^, the gra- 
dient of with respect to the parameters (Ao,mo) will 
be zero and the mass of ground state toq is estimated. 



B. Anisotropic Lattice Actions 

To examine the H-dibaryon in lattice QCD, we explore 
the improved actions on anisotropic lattices. These ac- 
tion display nearly perfect scaling, thus lattice-spacing 
artifact contributions are expected to be small, and pro- 
viding reliable continuum limit results at finite lattice 
spacings can be obtained. With most of the finite-lattice 
artifacts having been removed, one can use coarse lattices 
with fewer sites and much less computational effort. Us- 
ing a tadpole-improved anisotropic gauge action , we 
generate quenched configurations on a 12'^ x 60 lattice 
at five couplings in the range (3 = 2.0 — 4.0 and at a 
bare anisotropy of ^ = 5.0. We generated 500 gauge field 
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configurations for each lattice and the configurations are 
separated by 100 compound sweeps after skipping 1000 
sweeps for the thermahzation. We define a compound 
sweep as 5 over-relaxation [2^ sweeps followed by one 
Cabbibo-Marinari [2^ sweep. 

For the fermion fields, we employ the space-time asym- 
metric clover quark action on anisotropic lattice [23, 
with spatial Wilson parameter = 1. The clover im- 
provement coefficients Cs^t are estimated from tree-level 
tadpole improvement whereas for the ratio of hopping 
parameters C, = Kt/ Kg we adopt both the tree- level im- 
proved value and a non-perturbative one. Since it be- 
comes harder to obtain a reasonable signal-to-noise ratio 
at lighter quark masses for the multi-quark systems, we 
employ relatively heavy quark masses in our calculations. 
The bare strange quark mass is set by measuring the 
ss pscudoscalar mass at four heavy quark hopping pa- 
rameters Kh- At each strange quark mass, hadron prop- 
agators are measured for six light hopping parameters 
Ki such that the mass ratio of Mk/Mn compares well 
with the experimental value. Our quenched quark prop- 
agators cover a range of quark masses, corresponding to 
pion masses from 1325 MeV down to 500 MeV. We also 
considered two smaller masses, but find that the signal 
for these becomes highly unstable, hence do not include 
these in our analysis. 

C. Smearing technique 

To increase the overlap of the operators with the 
ground state, all of the hadronic correlators were calcu- 
lated using the method of smearing the interpolating op- 
erator, essentially making the hadronic operator spread 
around their central location in space. In this study, we 
use the gaussian smearing which is obtained by replacing 
the quark field q(x) by the smeared quark field qsmeari^) 
defined as [29j 

gsmcar(t,2^) ^ iV ^ exp{- ~ f ' }q{t,y), (7) 

y 

where is an appropriate normalization factor and p the 
smearing size parameter. This technique has numerical 
advantages since the smearing function separates into two 
factors one belonging to the quark and the other to the 
antiquark, thus will help to maximize the ground state 
contribution relative to the ones of the excited states. 
The problem is that the smeared operators are no longer 
gauge-invariant because the quark and the antiquark are 
spatially separated. We employed Coulomb gauge fixing 
to overcome this problem. 

D. Extrapolation to the physical quark mass and 
Continuum limits 

Chiral extrapolations of the H-dibaryon mass and 
binding energy to the physical point are important issues. 



In the exact SU(3) flavour symmetry, the non-interacting 
I{J^) = 0(0+) with strangeness s = —2 ground state is 
multiple degenerate, comprised of the states AA, A^S and 
SS with H-dibaryon as the ground state. A tightly bound 
H-dibaryon would indicate the chiral expansion of the 
form of that for single hadrons. The chiral extrapolation 
of single hadrons, such as the lowest-lying octet baryon 
masses, is an ongoing topic of discussion and motivates a 
deeper understanding of extrapolation form. The baryon 
chiral perturbation theory seems reluctant to reproduce 
LQCD results for the octet baryon masses, including the 
results for nucleon mass. Leinweber et al [s^] demon- 
strated that the chiral extrapolation method based upon 
finite-range regulator leads to extremely accurate value 
for the mass of physical nucleon with systematic errors 
of less than one percent. 

To addressed the challenges of SU(3) chiral perturba- 
tion theory to describe the baryon masses, Walker-Loud 
et al., detailed a comprehensive chiral extrapolation anal- 
ysis of the octet and decuplet baryon masses, using both 
the continuum SU(3) heavy baryon xPT^ as well as its 
mixed action generalization (sTI . [s^l ■ The results placed 
signature of linearity of nucleon mass in , providing a 
remarkable agreement with both the lattice data as well 
as the physical nucleon mass. This is in contrast with 
the expectations of chiral limit expansion of the gen- 
eral form MN{m^) = a + bm'^ + O(m^), where a and 
b are parameters determined from the lattice QCD data. 
Sharpe and Labrcnz [33| also found a more complicated 
and available form of chiral expansion of baryon masses 
in quenched approximation. The problem is that we have 
only several binding energy for each lattice spacing, us- 
ing the form made by Sharpe and Labrenz with so many 
coefficients may provides an unavailable and unreliable 
fit. Considering the form of binding energy is not known 
yet, we apply the general one to perform the chiral ex- 
trapolation and in fact we obtain a good result. 

To avoid the ambiguity in the chiral limit estimates, we 
extrapolate mass difference and AM = Mh — 2Ma and 
mass ratios AM/M\ using the simplest ansatz consistent 
with leading order chiral effective theory, 

f^a + l3x, (8) 

where x represents the pion mass squared and a and f3 are 
fit parameters. The quantities / and x are accompanied 
by statistical errors. We intend to find the combination 
of a and l3 which minimizes 

^ {f{a,P;x.)-{f,)f 

where i indexes different data points {x, /} and a is the 
statistical error of each quantity. The extrapolation is 
taken to physical point at fixed strange mass and M\ is 
taken as experimental input to make physical predictions. 

The continuum extrapolation for the chirally extrapo- 
lated quantities is another important issue in lattice cal- 
culations. The possible error that might effect the sim- 
ulation results comes from the scaling violation for our 
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action. Expecting that dominant part of scaling violation 
is largely eliminated by tadpole- improvement, we adopt 
an -linear extrapolation to the continuum limit, since 
the lattice-spacing artifacts in our calculations are ex- 
pected to scale as 0{al). Also, since the 0{a'^) effects 
largely cancel in forming the binding energy, we expect 
such contributions to be small. 



III. RESULTS AND DISCUSSION 

Typical examples of the effective mass plot at (ki, (3) = 
(0.3110,3.60) and (0.3115,4.00) arc shown in Fig. [H As 
can be seen, smearing improves the overlap with the H- 
dibaryon ground state resulting in an earlier plateau. 
Consequently the contributions of excited states were 
substantially reduced. We find clear signals up to larger 
time separations with insignificant statistical fluctuation 
domination. The fit range [tmimtmax] is determined by 
fixing tmax and finding a range of tmin where the ground 
state is stable against tmin- The statistical error analy- 
sis is performed by a single-elimination jackknife method 
and the goodness of the fit is gauged by the per degree 
of freedom, chosen according to criteria that x^/^df is 
preferably close to 1.0. The resulting effective masses 
of H-dibaryon and A states for other values of atm^ at 
/3 = 3.60 are tabulated in Table ID 



enough that the contributions of excited states are sup- 
pressed. Using the Levenberg-Marquardt algorithm, we 
indeed found such a range of t where linear term suffices 
in the data presented here. Fig. [5] shows the mass split- 
ting between the H-dibaryon and 2A threshold for the 
parameter combination (0.3110,0.3115) at as = 0.211 
and 0.188 fm, respectively. In the time interval where a 
single state dominates, the plateau region is reasonably 
consistent with that obtained for the effective mass of the 
H-dibaryon. The energy gap shows the negative value in 
the plateau region of 12 < t < 22 and seems more pro- 
nounced with H-mass smaller than two A's. The signal 
of mass difference is dominated by the large fluctuations 
in the H-dibaryon correlators beyond t ~ 22. 



FIG. 2: Effective mass difference between the H-dibaryon 
state and the 5- wave A -f A two-particle state at (k; , /3) = 
(0.3110,3.60) (left panel) and (k;,/?) = (0.3115,4.00) (right 
panel). 



FIG. 1: Effective masses of the I(J^) = 0(0+) stranglet 
(solid triangles) and lambda baryon (solid circles) at {fti, P) = 
(0.3110,3.60) (left panel) and {ki,/3) = (0.3115,4.00) (right 
panel) . 



TABLE I: Effective masses of lambda baryon and H-dibaryon 
on the 12'^ x 60 lattice at /? = 3.60 for various values of atm-^. 



atrriTr 


at Ma 


atMh 


0.38847(77) 


0.5202(24) 


1.0320(71) 


0.36661(78) 


0.4979(24) 


0.9865(71) 


0.34327(78) 


0.4746(24) 


0.9386(72) 


0.31759(77) 


0.4489(24) 


0.8868(72) 


0.28848(77) 


0.4199(24) 


0.8287(73) 


0.25438(77) 


0.3861(24) 


0.7606(73) 



In order to determine the energy difference AM = 
Mh — "iMiA precisely, we work in a regime where t is small 
enough that tAM << 1, and at the same time t is large 



The results on the other lattice spacings shows consis- 
tency in the behaviour of mass difference over the range 
of our pion mass range (see Tables [III HH] and IIV|) . 



TABLE IL Mass differences and mass ratios between the H- 
dibaryon and (A -I- A) two-particle state for various values of 
ajm^ at /3 = 2.00. 





ajAM 


AAf/A/A 


0.57298(32) 


-0.0351(45) 


-0.0396(51) 


0.55564(31) 


-0.0348(42) 


-0.0401(49) 


0.53723(30) 


-0.0319(44) 


-0.0375(51) 


0.51770(29) 


-0.0335(43) 


-0.0402(51) 


0.49682(29) 


-0.0331(42) 


-0.0407(52) 


0.47430(28) 


-0.0312(42) 


-0.0394(52) 



With all prerequisites available to measure the energy 
shift of H-dibaryon relative to the 2A threshold, we dis- 
play, in Fig. 131 the resulting mass differences extrap- 
olated to physical quark mass value using the ansatz in 
Eq. ([U . We note that slope of a linear fit in is slightly 
different at all lattice spacings. On the other hand, mass 
difference is almost constant and weakly dependent on 
quark mass. Nevertheless the results on all lattice spac- 
ings exhibit a negative value in the physical region. The 
negative mass difference observed in this region would 
imply an attractive interaction and hence a signature of 
H-dibaryon as a bound state. 
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FIG. 3: Effective energy shift between the H-dibaryon state 
and A + A two-particle state as a function of atiri^ squared. 
Solid circles, diamonds and triangles show the results at /3 = 
2.00, 2.80 and 3.20, respectively. 



FIG. 4: Plot of mass ratio AM/Ma as a function of ajm^ 
squared. Solid circles, diamonds and triangles show the re- 
sults at /3 = 3.20, 3.60 and 4.00, respectively. Dotted lines are 
the linear extrapolations to the chiral limit. 



TABLE III: The same as Table [IT] but at /3 = 3.20. 



atm^ 


at AM 


am/Ma 


0.41147(45) 


-0.0125(49) 


-0.0223(87) 


0.39017(46) 


-0.0138(49) 


-0.0256(92) 


0.36742(46) 


-0.0140(49) 


-0.0271(94) 


0.34252(46) 


-0.0143(48) 


-0.0291(98) 


0.31463(46) 


-0.0141(49) 


-0.0304(104) 


0.28250(47) 


-0.0143(49) 


-0.0332(112) 



TABLE IV: The 


same as Table HIl but at /3 = 4.00. 


atirin 


at AM 


am/Ma 


0.37144(91) 


-0.0070(70) 


-0.0141(141) 


0.34893(91) 


-0.0094(70) 


-0.0199(148) 


0.32486(90) 


-0.0095(70) 


-0.0211(155) 


0.29802(89) 


-0.0098(70) 


-0.0233(165) 


0.26725(88) 


-0.0101(70) 


-0.0260(178) 


0.22738(111) 


-0.0109(71) 


-0.0307(199) 



Since the quenched spectroscopy is quite reliable for 
mass ratio of stable particles, it is physically even more 
motivating to extrapolate mass ratio instead of mass. 
This allows for the cancellation of systematic errors 
since the hadron states are generated from the same 
gauge configuration and hence systematic errors are cor- 
related. Fig. |4] shows the chiral extrapolation of the 
ratio AM/M\ at our smaller lattice spacings. The ratio 
shows a weaker dependence on quark mass and moves 
into the physical region with a trend of attractive inter- 
action. In the chiral limit the estimated mass difference 
at our smallest and largest lattice spacings are consistent 
with those obtained by NPLQCD and HALQCD Collab- 
orations [l^, [l^ . The uncertainties from the chiral fits 
range from 2 to 10 percent for to the smallest to largest 
lattice spacings explored here. 

The quenching effects might be one of the largest 
source of the systematic uncertainties. However, with 
appropriate definition of scale, the mass ratios of stable 
hadrons differ from the corresponding observed values by 
less than 6% in quenched approximation [13, Ull . In or- 
der to absorb as many quenching effects as possible, we 
set the scale by the physical Kg by calculating the ra- 
tio Ma/Mn. We found that the ratio deviates in the 



range of 3 — 4% from its experimental value verifying 
that the value of Kg used is very close to the physical 
quark mass. However, since our calculations use rather 
heavy pion masses (> 500 MeV), the quenching effects 
are less noticeable within our statistics. We include a 
modest estimate of order 5% quenching uncertainties in 
our analysis. 

Whether the energy difference moves in the continuum 
with an attractive interaction needs to be explored. Since 
the finite-spacing errors in our calculations are expected 
to scale as a^, we expect such contribution to have a 
small effect on binding energy. Consequently, we expect 
the observation of H-dibaryon to survive the continuum 
extrapolation. We perform the continuum extrapolation 
of the chirally extrapolated mass ratios in Fig. [S] and 
present the results in Table fVl Using an a^-linear extrap- 
olation, we adopt the choice which shows the smoothest 
scaling behaviour for the final values, and use other to 
estimate the systematic errors. 

As is clear from the figure, the mass ratio shows a 
weak dependence on the lattice spacing and varies only 
slightly over the fitting range. Thus we expect our con- 
tinuum extrapolation accurate and unambiguous. The 
continuum extrapolation is accompanied with an order 
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FIG. 5: Compilation of results for the mass ratio, AM/Ma 
in the continuum limit. The solid line represents linear 
extrapolation to the physical limit. 

TABLE V: Mass ratios between the H-dibaryon and (A + A) 
two-particle state at various lattice spacings. 

a^(fm) am/Ma 

0.3974(34) -0.0409(66) 
0.3234(43) -0.0391(73) 
0.2599(62) -0.0486(78) 
0.2347(57) -0.0422(149) 
0.2114(70) -0.0403(242) 
0.1875(56) -0.0392(296) 



possibility of the volume dependence of the binding en- 
ergy in the large volume limit. We intend to pin down 
this problem for a conclusive signature in future work. 



IV. CONCLUSIONS 

The question of whether the H-dibaryon is a bound or 
unbound state is still under debate. The observed nega- 
tive energy shift pattern favours the former. In conclu- 
sion, we have presented evidence for existence of bound 
H-dibaryon state in the physical limit from quenched lat- 
tice QCD calculations. The calculations were performed 
over a range of pion masses and lattice spacings using 
improved anisotropic lattices with refined computational 
techniques. Attractive interaction was found in chiral 
limit for all pion masses used in this study and the con- 
tinuum limit estimate seem to agree with the predicted 
value, which is one of the main results of our paper. Our 
results seem to be consistent with the recent results of 
NPLQCD and HALQCD Collaborations for the energy 
shift. Our analysis takes into account possible artifacts 
such as, statistical, chiral and continuum extrapolation 
uncertainties and those arising from quenching effects. 
On the basis of our lattice calculation we speculate that 
the H-dibaryon is to be identified as bound state. How- 
ever, the final conclusions will have to await dynamical 
simulations incorporating a systematic study of various 
possible interpolators that are likely to have a good over- 
lap with H-dibaryon. We plan to further develop this 
calculation to involve a combination of AA — SS — EE 
interpolators on larger volumes. 



8% systematic uncertainty from the linear fit in a^. Us- 
ing the physical A mass M\ = 1115.68 MeV, we obtain 
a continuum estimate of binding energy 47 ± 37 MeV 
for the binding energy. The uncertainty shown here re- 
sults from the statistic and systematic uncertainties com- 
bined in quadrature. The systematic uncertainty includ- 
ing quenching, chiral and continuum extrapolation effects 
is estimated to be of the order of 15%. Note that we can- 
not estimate the finite size effects since we have been 
working with one lattice size. Even though our spatial 
extent of L is reasonably large, we cannot rule out the 
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